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1. Introduction 

Ultra-high energy neutrinos and their interactions continue to attract much atten- 
tion. This is in spite of the fact that no such neutrino has been seen so far (and 
hence bounds been placed on their fluxes [1]). Much of the continuing interest has 
been occasioned by the observation, in more than one detector, of ultra-high energy 
cosmic rays (UHECR). An interaction of such cosmic rays with either the microwave 
background radiation or even the atmosphere would presumably lead to the genera- 
tion of charged pions and through their decay, to extremely energetic neutrinos [2]. 
Alternatively, primary ultra-high energy neutrinos themselves could lead to UHECR, 
thereby avoiding the GZK bound. A possible source for such primary neutrinos is the 
decay of an extremely massive primordial relic or even a cosmic string [3] . Whatever 
their origin, it can safely be asserted that cosmic high energy neutrinos are inex- 
tricably linked to the very high energy cosmic rays. The experimental detection of 
such neutrino fluxes is thus expected to provide rare insight into the origin of such 
cosmic rays and probably to physics beyond the SM as well. Such observations have 
the additional promise of probing stellar structures [4] , for unlike charged particles, 
cosmic neutrinos do not suffer any bending due to inter-galactic magnetic fields and 
hence arrive on earth in a direct line from their source. 
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Consequently, various experimental efforts are being planned. Pilot experiments, 
based on the optical detection of Cerenkov light emitted by the muons created in 
charged current reactions of neutrinos with nucleons either in water or in ice, include 
the Antarctic Muon And Neutrino Detector Array (AMANDA) [5] in the South Pole 
ice and the one at Lake Baikal [6]. The next generation experiments using simi- 
lar techniques comprise the Neutrino Telescope and Abyss environmental RESearch 
(ANTARES) [7], the Neutrino Experiment SouthwesT Of GReece (NESTOR) project 
in the Mediterranean [8], as well as ICECUBE [9], the proposed kilometer scale ver- 
sion of the AMANDA detector. Recently, arguments have been forwarded in favour 
of facilities based on the detection of radio pulses emanating from the electromagnetic 
showers created by neutrino interactions in ice and other materials. The primary ad- 
vantage of such a technique would be the scalability up to an effective area of 10 4 km 2 
and the Radio Ice Cerenkov Experiment (RICE) experiment at the South Pole [10] is 
a functioning prototype. It has also been realized that neutrinos can initiate horizon- 
tal Extensive Air Showers (EAS) which could be detected by giant ground arrays and 
fluorescence detectors such as the cosmic ray Pierre Auger Project [11]. Deeply pen- 
etrating EAS could also be detected by observing their fluorescence light from space 
based instruments such as the Orbiting Wide-angle Light-collector (OWL) [12] and 
the Extreme Universe Space Observatory (EUSO) [13]. Finally, there is the newly 
approved balloon experiment ANtartic Impulsive Transient Antenna (ANITA) which 
will look for radio Cerenkov pulses created by ultra-high energy neutrino interactions 
and emanating from very large chunks of the Antarctic ice cap. Its energy threshold 
is about 10 18 eV, so it will primarily be looking for GZK neutrinos [14]. 

These experiments, taken together, are sensitive to neutrino energies of upto 
10 11 GeV or so. The actual event rates are somewhat uncertain though, as they 
depend crucially on both the predicted neutrino fluxes, as well as on the ultra- 
high energy neutrino cross sections which, for want of a better method, we may only 
estimate by a reasonable extrapolation beyond the measured regime. The interaction 
of UHE neutrinos with matter is through deep inelastic scattering of the neutrinos 
with protons and neutrons. Over the last few years, numerous issues with regard to 
the nature of the cross section of v — N scattering (where N is a proton or a neutron) 
have gained importance and some of these are discussed in Refs. [15-17]. 

Most of the discussion on v — N scattering has been based on the leading order 
(LO) expressions for neutrino nucleon scattering (i.e. a s independent). The usual 
procedure followed has been to use the lowest order parton level cross section and 
convolute it with the LO or sometimes even next to leading order (NLO) parton 
distributions. QCD corrections to the partonic cross sections have typically been 
neglected, in view of the high energies involved and the consequent small value 
of the strong coupling constant a s . While, at first sight, such an approximation 
may seem appropriate, it must be borne in mind that the consequent uncertainties 
may limit the sensitivity of neutrino telescopes [5, 7-10] (and to a smaller extent, 
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the cosmic-ray detectors [11-13]) to physics beyond the Standard Model. Amongst 
possible such scenarios, of particular interest are theories with supersymmetry [18], 
extended gauge or higgs sector [19] or, more recently, those with a low energy gravity 
sector [20]. Perhaps, of even more importance, are the effects on the determination of 
neutrino mixing parameters [21], and neutrino-tomography of the earth's interior [22]. 
The importance of a more accurate estimation of the neutrino interaction rates as 
well as their kinematical distributions, thus, cannot be overstated. 

In this paper, we explicitly calculate the 0(a s ) QCD corrections to the parton 
model result and show that while it is not very large, it is by no means negligible. 
Moreover, we study carefully the behavior of this correction as a function of neutrino 
energy and find behavior which is not necessarily very intuitive. For example, there 
is a delicate interplay between the magnitude of a s , the structure of the higher order 
integrals and the size of the parton distributions (particularly the gluon) in LO and 
NLO. This gives a non-trivial energy dependence to the ratio of the LO and NLO 
cross section (which we will call the if-factor). We work throughout in the MS 
scheme. 

Another issue which is of relevance in these energy ranges (and which, again, has 
been addressed in Refs. [15-17]) is the question of carrying out perturbative calcula- 
tions at ultra low Bjorken x (down to 10~ 8 ). No data exists to help in parametrisa- 
tions of parton distribution functions at such values of x and one can only be guided 
in these regions by a somewhat improperly understood physical picture of a highly 
dense nucleon of partons. Data from HERA stops around x ~ 10~ 5 or so and below 
that, some physical picture of shadowing and saturation effects (particularly at low 
Q 2 ) needs to be incorporated to understand the physics of a nucleon with a high 
density of partons. 

We have tried to address both these issues in this paper. We have explicitly 
calculated the 0(a s ) corrections to the partonic cross sections and convoluted them 
with appropriate parton distributions. We have also addressed the issue of extrapo- 
lation of the partonic distributions to regions where simple DGLAP evolution is not 
expected to hold. 

The paper is organised as follows. In Section 2, we present a discussion and 
justification for the various partonic distributions that we have used in various parts 
of the (Q 2 ,x) plane. In Section 3, we present detailed expressions for the 0(a s ) 
corrections to the lowest order partonic cross section for neutrino and antineutri- 
nos scattering against an isoscalar target. These expression are, of course available 
elsewhere but for the sake of completeness and clarity, we feel it would be useful to 
present them in a form that is amenable to discussions later in this paper. In Section 
4, we present our results for LO and NLO cross sections, both for the differential 
distributions da/dlogx and da / d log Q 2 , as well as the total cross section. This has 
been done for neutral as well as charged current cross sections. Section 5 has a short 
discussion on saturation and the final section makes a few concluding remarks. 
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2. Parton Distributions 

As already mentioned briefly in the introduction, deep inelastic scattering of UHE 
neutrinos is unique in the sense that it explores extreme regions of the (Q 2 , x) phase 
space where no data from terrestrial accelerators exist to help with the nature of the 
distributions. We are therefore forced into making some assumptions regarding the 
nature of these distributions, particularly in the region of ultra low x and low Q 2 
(the so-called saturation region), as well in the region of ultra low x and very high 
Q 2 , for which too, no data exist. In what follows we shall expand a bit on this theme. 

In the region between, say 50 GeV 2 < Q 2 < 10 5 GeV 2 , standard parametrisations 
are expected to work reasonably well. All parametrisations in these regions use 
parton densities of the form Ax a (l — x) /3 f(y/x) where the last function is some 
polynomial in y/x. We have checked explicitly that all three of CTEQ5 [23], MRS99 
[24] and GRV98 [25] parametrisations give very similar results within this region. 
However, the explicit parametrization for either of CTEQ and MRS ditributions 
do not work below x < 10~ 5 whereas GRV98 allows x values upto 10~ 9 . For our 
numerical results, we therefore work with GRV98, except in the saturation region 
and the very large Q 2 and small x region. 

The region above Q 2 = 10 6 GeV 2 and low x (x < 10~ 5 ) is beyond the re- 
gion of validity of the explicit parametrisations of GRV (which is valid only upto 
Q 2 = 10 6 GeV 2 ). One option here is to use the starting distributions for GRV and 
evolve them (to LO or NLO order) into the above region. However, in this region, 
one also has the option of using an analytical form for the distribution function - the 
so-called double leading log (DLL) or Double Asymptotic Scaling (DAS) forms. In 
this approach, the DGLAP equations can be solved analytically, assuming that the 
evolution is driven by the splitting function P gg . This was first shown in [26] and 
later developed for ultra-high energy neutrinos in [27]. A general review with calcula- 
tional details may be found in [28]. This region is then matched appropriately at the 
boundary with GRV. One could also use, other than variations of the DGLAP ap- 
proach, the ln(l/x) resummation programme of BFKL or the unified DGLAP-BFKL 
treatment of [37]. The numerical changes are insignificant. 

For the other region where the explicit GRV parametrization is not valid, namely 
large Q 2 (> 10 6 GeV 2 ) and moderate x (x > 10~ 5 ), one may again evolve the parton 
densities appropriately and then use them. However, numerically, it is simpler to 
just use GRV98 with a fixed value of Q 2 = 10 6 GeV 2 , and, within this region, the 
error due to such an approach is negligible. It might seem counterintuitive that the 
final result depends negligibly on which resummation scheme is used in the ultra-low 
x and/or very high Q 2 scheme. However, as we show later in this paper, this simply 
follows from the fact that at these extreme (x, Q 2 ) values, the relative contribution 
in the integral (over x and Q 2 ) to the total cross section is already very small. 

The summary of the various evolution equations and their applicable regions 
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Figure 1: Summary of the parton distributions used in various regions DAS stands for 
double asymptotic scaling (or equivalently DLL or double leading log) 

are shown in Fig. |l| We have used the leading log expressions for DAS; this is 
appropriate for matching with NLO GRV. 

Some comments regarding the saturation region are in order. In DIS at very low 
x and low Q 2 , we find a high density of partons which is a non perturbative system 
even though a s may be small. There are essentially two approaches to high density 
QCD - the GLR [29] and Mueller, Qiu [30] aproach and the effective Lagrangian 
method of McLerran and Venugopalan [31]. Using these two approaches a non-linear 
evolution equation has been developed by the Tel- Aviv group and a review may be 
found in their paper [32]. We have used saturated densities in the low x and low 
Q 2 region using this approach [33] . Some comments on the saturation region can be 
found later, in Section 5. 



3. NLO expressions for v — N scattering 

In this section we briefly review the 0(a s ) (NLO) expressions for neutrino nucleon 
scattering. Many of these details may be found scattered in the literature — for 
example, see [39] — however it is convenient to write them down here in some detail, 
clearly demarcating the LO and NLO pieces. We will assume an isoscalar target N 
[= (p + n)/2] as this is an excellent approximation to the typical neutrino detector 
material. 
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The general expression for the neutral current (NC) reaction 
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where Gf is the Fermi coupling, Mz the Z boson mass and s the total centre-of-mass 
energy squared. The quantity Q 2 (= — q 2 ) is the momentum transfer, the Bjorken 
variable x = Q 2 /(2p.q) with p (p 2 = M 2 ) being the four momentum of the proton 
(or the isoscalar target) and y = Q 2 /(xs). 

The structure functions Fj = F^x, Q 2 ) (i = 1, 2, 3) are given by the convolution 
of bare parton distributions qf\y p ) and the partonic structure function Fj. In other 
words, 

F i (x,g 2 ) = g/°)®Ff + 
where <S> denotes the convolution defined as 

" J dyp 

Denoting the generic weak interaction vertex involving a quark / by 7 M (cy ± 0^75), 
the leading order (LO) partonic structure functions are given by 

F t ?(z)=K i 6(l-z) (3.4) 

where 2Ki — K 2 — [{cy) 2 + (c^) 2 j and K 3 = 2 c v c* A . Clearly, K 1>2 correspond to 
parity conserving interactions whereas K 3 is a measure of parity violation. In the 
case of DIS, NLO corrections implies corrections to the parton structure functions 
coming from 0(a s ) corrections to the partonic cross sections. 

For the quark initiated process, to order a s , the partonic structure functions are 
given by 
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where the first term is the zeroth-order piece and the rest comprise the next-to- 
leading order corrections. The splitting function P qq is defined through 



PM = 



1 + z 2 



(3.6) 



where the "plus prescription" is defined as usual. The renormalisation scale \i is 
introduced to make the coupling dimensionless in n(— 4 + e) dimensions and 
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Clearly, the collinear singularity manifests itself in the ( dependent term in eqn. (|3.5|) 
in the limit e — > 0. And, finally, the coefficient functions Cf(z) are given by 
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A\ = -2z A q 2 = A q 3 = -{l + z), 

For the gluon initiated process, on the other hand, the parton structure functions 
are given by 

Q 2 - 



F?(z, Q 2 ) =K^ 
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Af = 4z(l -2) Af = A3 = 0. 



(3.10) 

(3.11) 
(3.12) 
(3.13) 



Note that the gluonic contribution to F 3 vanishes identically. This feature follows 
from the V — A structure of the weak interaction whereby the quark contribution 
is exactly opposite to the anti-quark contribution. For the anti-quark initiated pro- 
cesses all the structure functions are the same except for F q (z, Q 2 ) which appears 



with the opposite sign in eqn. (|3.1| ). Absorbing the mass singularity into the bare 
parton distributions in the MS scheme, the Q 2 dependent structure functions can 
be expressed in terms of the renormalised parton distributions through 



Fi(x,Q 2 



3=1 



+ a s {Q 2 



Cf(z) ® ( qj + qj) + 2Cf ® g 



, (3.14) 
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, (3.15) 
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F 3 (x,Q 2 ) = f>M {(^*,Q 2 )-^Q 2 )) 



+ a s {Q 2 



C&z) ® (g, - 



(3.16) 



In eqns.( 3.14 - 3.16| ), the first term refers to the LO contribution to Fj(x, Q 2 ) whereas 
the second term is the pure NLO contribution. The superscripts pc and pv stand for 
'parity conserving' and 'parity violating' respectively. The fact that the tree-level 
relation F 2 = 2xFi is no longer valid will play a key role in understanding some of 
the results detailed in the next section. 

It is useful, at this stage, to note that for an isoscalar target we have 



1 



+ K { / c) (s + s + b + l) +K^(c + c + t + tj, (3.17) 
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where K u and Kd are the couplings of the up and down type quark respectively. All 
parton distributions are the "renormalized" parton distributions (i.e. q, = qj(x,Q 2 )). 
The general expression for the charged current (CC) reaction 
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(3.18) 



all the notation being the same as NC. As the couplings Ki are now flavour indepen- 
dent (2Ki — K-2 — K 3 = |), they have been absorbed. For the isoscalar target, the 
quark distributions have the form 



E? (<!, '/, 



it + d + u + d 
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4. Results and discussion 



Having set up the formalism, we now turn to an estimate of the numerical size of 
the NLO corrections. While the exact values of the cross sections for neutrino and 



- 8 - 



antineutrino induced processes differ, (as also for charged-current and neutral-current 
interactions), one expects some similarities in the relative size of the corrections. 

The total v and v cross sections are displayed in Fig. The rationale for choosing 
the particular scaling will become apparent as we progress. While the behaviour is 
essentially the same as that obtained by previous authors, the exact numbers (for the 
LO) differ marginally. The change is the result of our using a different set of parton 
distributions as compared to, say, Ref. [16]. Moreover, we have taken care to use the 
parton distributions appropriate for the order to which the partonic subprocesses are 
being calculated in each instance. To obtain a quantitative estimate of the differences 
in the total cross sections, it is useful to express them as a function of the neutrino 
energy. While a polynomial fit for log a in terms of log-E is straightforward, it is, 
unfortunately, not very illuminating. We attempt, instead, a piecewise fitting of the 
form 

with the parameters A and 7 as given in Table |l|. Thus, while a grows almost linearly 
with E u for relatively small neutrino energies, the growth is tempered to ~ E® A for 
larger E u . The particular scaling in Fig. |^, thus, serves to highlight the deviations 
at high energies from the exact scaling relation of eqn. (|4.1|) . Such a behaviour for 
the total cross section was found earlier, in [34], and as demonstrated there, it is 
the result of the fact that when the differential cross section is integrated over x 
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Figure 2: The total NC and CC cross sections as a function of the (anti-) neutrino energy. 
The left and right panels correspond to neutrino and antineutrino scattering respectively. 
The solid (dashed) line represents the NLO (LO) cross sections. 



and Q 2 , the Q 2 integration is effectively cut-off at Q 2 ~ M\ w coupled with the 
fact that the number of valence quarks is finite. One might wonder at this stage 
at the apparent lack of unitarity reflected by the value of the exponent 7 at high 
energies. In particular, in view of the rapidly rising parton distribution functions 
at low x predicted by all evolution equations, the question of the total cross section 
saturating the unitarity bound already at around E v ~ 10 8 GeV is a serious one that 
is addressed, for example, in Refs. [15,17]. 
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NLO 


2.71 x 1CT 3 
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0.360 


D-NC 


LO 


1.44 x 1CT 3 
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1.16 


0.391 


NLO 


1.52 x 1CT 3 


0.963 


2.06 


0.360 


v-CC 


LO 


8.89 x 1CT 3 


0.940 


3.57 


0.383 


NLO 


8.82 x 1CT 3 


0.941 


6.50 


0.349 


D-CC 


LO 


3.95 x 1CT 3 


0.961 


2.85 


0.386 


NLO 


4.15 x 1CT 3 


0.960 


5.18 


0.353 



Table 1: Parameters for the piecewise single power fit (see eqn.(4-l)) for the cross sections 
in terms of the ( anti- Jneutrino energy. 

While both Fig. and Table [I] seemingly demonstrate that, for small neutrino 
energies, the difference between the LO and the NLO cross sections is miniscule, the 
same certainly cannot be said of the results at large values of E u , This difference can 
be quantified in terms of the i^-factor: 

K = (4.2) 

which we plot in Fig. ^. 

Various points are worthy of note here. As already anticipated, for small neu- 
trino energies, the i^-factor is remarkably close to one and confirms the usual belief 
that QCD corrections are unimportant. At larger energies [E v > 10 4 GeV), the 
NLO contributions start to become important and the .ff -factor slowly grows to its 
maximum (~ 1.08) at about E v ~ 2 x 10 6 GeV. Thereafter, the onlo grows slower 
with energy than does cxlo (see Table HD, resulting in a rather strong reduction of 
the i^-factor by the time one reaches E v ~ 10 10 GeV. 
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Figure 3: The ratio of the NLO to LO total NC and CC cross section as a function of 
( anti- jneutrino energy. 

4.1 z/-NC scattering 

The reasons for the afore-mentioned behaviour of the i^-factor are manifold and 
intertwined. Hence, before we attempt an understanding of this, it is perhaps more 
useful to look at the differential distributions. We begin by concentrating on the 
neutral current scattering process of neutrinos off isoscalar targets. In Fig. [|, we 
show the variation with Q 2 for three representative neutrino energies. The shape of 
the distributions are governed by both the dynamics as well as the parton densities. 
For low neutrino energies, the maximal energy transfer falls well short of Mz and thus 
the neutrino-parton interaction is well described in terms of a 4-fermi interaction. 
As is well known, the cross section due to such an interaction Lagrangian grows with 
the available center of mass energy. 

On the other hand, the partons densities fall very sharply as the fraction of the 
proton energy they carry approaches unity. The interplay of these two effects leads 
to the skewed bell-like shape of the distribution. As the neutrino energy increases, 
so does the typical value of the neutrino-parton center of mass energy. Once this 
crosses s ~ (lOOGeV) 2 , the natural scale of the problem, the contribution to the cross 
section is naturally dominated by virtual exchanges with Q 2 ~ M|. The skewness 
of the distribution is therefore removed to a large extent (see the central panel of 
Fig. At even larger energies, this tendency is only reinforced. 

The issues discussed above can also be understood in terms of distributions in 
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Figure 4: The differential NC cross section wrt log(Q 2 ) for three different neutrino ener- 
gies. In each case, the solid (dashed) line give the NLO(LO) cross sections. The quantity 
Q gives the factor by which the graphs have been scaled. 



the Bjorken variable x (see Fig. |^). With the fast fall of the parton densities at large 
x, it is only natural that such regions of the phase space should contribute very little 
to the cross section. This, of course, is reflected in the first two panels of Fig. || 
However, for very large E u , the situation is changed somewhat. Dynamics dictates 
that the region Q 2 ~ M§ (and, hence, relatively small x) receives prominence in 
the integration. As the maximum in distribution is still far away from x < 1, the 
aforementioned damping does not come into play and the fall-off is much smoother. 

With these features in mind, we now return to the issue of the i^-factor. To 
start with, let us re-examine Fig. |], more particularly the last two panels where the 
deviation is more pronounced. Even a casual glance (and this is indeed borne out by 
a quantitative test) shows that the ratio of the NLO and LO differential cross sections 
is roughly constant over the significant range of Q 2 . The shift from unity is as large as 
25% for extremely large E v . As for the x distribution, a similar effect is in operation 
for not too large neutrino energies (E u < 10 8 GeV). For progressively larger E^s 
though, the effect is concentrated more towards the peak in dcx/dlogx (and may be 
as large as ~ 30%). We may thus conclude that, while the NLO correction tends to be 
nearly independent of Q 2 , it certainly does have a non- negligible variation in x, with 
the largest shifts occurring for x < M§ / (2M P E U ) . In other words, were we to consider 
the normalized differential distributions, the LO and NLO curves for da/dlogQ 2 , 
for a given E v , would be virtually indistinguishable from each other, while those for 
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Figure 5: As in Fig.fy but for variation in log(x) instead. 



dcx/dloga; would seem laterally displaced, with the effect more pronounced for larger 
E v . 

The various contributions to the If -factor are often separated (albeit not in a 
strictly gauge-invariant way) into two parts: those arising from the differences in 
the LO and the NLO parton distributions and those that are due to the changes in 
the partonic cross sections. In Fig. ^ we display the effects of such a division. This 
would also serve to facilitate a comparison of our results with those available in the 
literature. As is seen readily, ignoring the corrections to the lowest order partonic 
cross sections but using the NLO distributions would have led to a much steeper 
suppression for larger neutrino energies. The main reason for the fall is that while 
at large E v and hence large s, x is small and hence the partonic distributions are 
large, the rise of the distributions is steeper for LO than for NLO. (This is, of course, 
a well known fact and arises because, at NLO, the splitting function (say P qg ) has a 
1/x singularity at 0(a s ) which needs to be stemmed by choosing flatter distributions 
at a given energy compared to the LO case. Put differently, a stable evolution of, 
say, F 2 requires that the NLO steepness of P qg has to be compensated by a gluon 
density which is less steep at NLO than at LO. It is for this reason that the starting 
distributions for LO and NLO are always substantially different). In the literature 
(for example, Ref. [16]), it has often been claimed that using the LO partonic cross 
section convoluted with the partonic distributions in the DIS scheme is an effective 
method for approximating the magnitude of the full NLO corrections. However, as 
Fig. 6 clearly shows this is not the case and the 'true' NLO corrections (wherein the 



NLO partonic cross section is convoluted with the NLO distributions), gives a result 
very different from using just the LO partonic cross section convoluted with either 
the NLO MS distribution or the distributions in the DIS scheme. 
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Figure 6: A comparison of the true K-factor (solid line) for the neutrino-nucleon neutral 
current process and that obtained if the leading order cross sections were convoluted with 
the NLO (MS) parton distributions (dashed line) or the DIS ones (dotted line). 

Having understood a decreasing i^-factor for high neutrino eneries, we now turn 
to the presence of a maximum. For sufficiently small E v 's, when x is not too small 
(see Fig. the gluon, which is what may make the NLO large at small x, is still 
to kick in, and the LO and NLO curves for the total cross section start out virtually 
together. As we move to larger values of E v , and hence smaller x, the gluon density 
itself becomes larger. Consequently, onlo increases faster than ctlo- F° r even higher 
Ey's, the bulk of the contribution comes from progressively smaller x and the fact 
of the LO distributions being larger, at small x, than the NLO ones becomes the 
overriding factor; the rise in ctnlo is suppressed. 

Before concluding this section, we would like to recall that whereas the gluon 
is the main contributing factor in the NLO graphs, no gluon appears in the initial 
state in the LO case. However, even for the LO, the gluon manifests itself, indirectly, 
through large values of the sea quark distributions. The inclusion of the NLO con- 
tribution only serves to ameliorate this suppression. It is interesting to note that 
the 'extra' piece may contribute as much as 10% to the cross section and is, by no 
means, negligible. 
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4.2 The z/-NC and the CC cross-sections 

Having discussed the z/-NC interaction in detail, we now comment on the other 
interactions that we have briefly discussed above. As both Fig. ^ and Fig. |3| show, 
the main features are quite similar and it is the small differences that we concentrate 
on here. Starting with the z/-NC process, note that the i^-factor at large energies is 
nearly the same as that for the corresponding z/-NC interaction. While the value at 
the peak (E D ~ 10 6 GeV) is also approximately the same, the dip at lower energies 
is far more gradual. Even more striking, at low energies E D < 10 GeV) is that the 
.fT-factor is discernibly different from unity. This difference is but an example of 
the afore-mentioned subtle interplay between the parton-level cross-sections and the 
parton densities. At the partonic level, the only difference between the z/-NC and 
the I/-NC cross sections lies in the reversal of the sign of the structure function F 3 , 
which in turn, is controlled by the valence quarks. At small energies, most of the 
contribution to the total cross section comes from the valence quarks and hence this 
extra additive piece has a significant role to play (and, indeed, is responsible for the 
difference in z/-NC and I/-NC cross sections). This relative shift is less pronounced 
for the NLO, as the F 3 contribution is now smaller while the combined Fi and F 2 
contribution is marginally larger - again, a consequence of the interplay between the 
parton-level cross sections and the shape of the parton densities themselves. With the 
increase of the (anti-)neutrino energy the peak in da/dlogx shifts to progressively 
smaller values of x (see Fig. [5[), thereby rendering the F3 contribution insignificant. 

As for the charged current processes, it is easy to see that the NLO corrections 
are quite analogous to the cases already considered. Thus, notwithstanding the 
difference in the absolute value of the cross sections, the i^-factor is expected to be 
similar, and this indeed is attested to by Fig. [| There is a minor difference though. 
Unlike the case for the NC processes, at high energies, the i^-factors for the neutrino- 
and antineutrino-induced CC processes are no longer the same. This has its origin 
in the expression for the F 3 term as applied to CC (see eq. ( p.!9| )). At high energies 
(and, hence, low x), the F 3 contribution is now nearly proportional to the strange 
sea content in the nucleon. As the latter is, by no means, small at such low x values, 
the differences in the absolute values of the cross sections as well as in the i^-factors 
are expected. 

4.3 Observables 'stable' under QCD corrections 

We have seen that all the four relevant cross sections receive significant corrections at 
the NLO level. Moreover, the corrections are not easy to parametrize as a function 
of energy. Added to this are dependencies on the parton distributions used. It is 
therefore of interest to inquire if there exist observables that are relatively insensitive 
to such corrections, not the least for their usefulness in searches for new physics. One 
such observable is the ratio of a^ c and a? c (the corresponding ratio for the NC cross 
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sections is not a measurable for cosmic neutrinos). In Fig. [7]a, we exhibit this ratio as 
a function of the (anti-)neutrino energy. Since the bulk of the NLO corrections are 
common to the two cross sections under discussion, the ratio, as expected, changes 
little when these corrections are included. 




10 io 3 io 5 io 7 io 9 io" 10 10 3 io 5 io 7 io 9 io" 

E v (GeV) E v (GeV) 

Figure 7: (a) The ratio of the charged- current cross sections for a neutrino and an an- 
tineutrino scattering off an isoscalar nucleon; (b) the ratio of charged- current and neutral 
current cross sections for a neutrino, an antineutrino and a beam composed of equal num- 
ber of neutrinos and antineutrino s. In either case, the solid (dashed) lines correspond to 
the NLO (LO) cross sections. 



Unfortunately, in the context of an actual experiment, the ratio above would 
be measurable only in the presence of a substantial magnetic field, and is hence of 
limited use as far as a neutrino telescope is concerned. However, the ratio of charged 
current and neutral current cross sections is measurable, and can indeed be a good 
discriminator in the search for new physics [18]. As Fig. [7]& shows, this ratio, again, 
is far less sensitive to QCD corrections than are the individual cross sections. 



5. Saturation effects 

We have seen in the earlier section the complex interplay between partonic cross 
sections and the evolution of partonic densities at LO and NLO. 

However, there is another issue which is of importance in this process. At the 
extremely low values of x that we are considering, the proton is a dense collection 
of partons and saturation effects are expected to come into play. Typically, this 
happens when recombination effects from g g — > g start having a discernible effect. 
Such processes are expected to slow down the usual rise of the gluon distribution 
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predicted by DGLAP evolution and is expected to take place around a scale 



The parameters A and xq as fitted to HERA data are found to be A = 0.288 and 
x = 3.04 x 10~ 4 [17,35]. Around x ~ 10~ 8 this works out to Q 2 ~ 20 GeV 2 . As is 
clear from the above, at larger values of x, saturation is expected to set in at even 
lower scales. 

Keeping this in mind, we have smoothly interpolated, in the low Q 2 ( < 50 GeV 2 ) 
and low x (< 10 -6 ) region, from GRV to the saturation model mentioned in the 
introduction [32]. All the graphs and the analyses discussed in this paper are based 
on interpolation of densities shown in Fig. 1. 

While this is, of course, the consistent way of doing this, a useful exercise would 
be to estimate the effect of these saturation models on our results. To this end, we 
substituted the saturation model in the low x and Q 2 region by the standard GRV 
evolution and compared the results. We find that the effect of saturation, in the 
total cross section is negligible and, at its largest, is of the order of around 0.2%. 
The reason for this is not far to seek. In the large range of x and Q 2 values that we 
are spanning in our integrals, the saturation region accounts for an exceptionally tiny 
part of the phase space and consequently, has little effect on the final result. The 
physical requirement of saturation is however, an important one for many reasons, not 
least in stemming the rise of the distributions at ultra low x and thereby preventing 
the cross section from crossing the unitarity bound. As mentioned earlier, we will 
not discuss this aspect here, but refer the reader to Refs. [15, 17]. 

6. Conclusions 

In this paper, we have explicitly calculated the NLO (i.e 0(a s )) corrections to ultra 
high energy neutrino-proton (or isoscalar nucleon) scattering. We have identified the 
various regions in the x-Q 2 space where different parametrisations of the quark and 
gluon densities are valid and discussed the contributions of each of these regions. 
In particular, we have carefully looked at the low Q 2 and ultra low x region and 
discussed the effect of saturation in this region and found it to have a very small 
effect on the overall NLO behaviour. We have also calculated the size of the NLO 
corrections compared with the LO (the K factor) and tried to give a qualitative 
explanation for the results. 

It is clear from our analysis that the higher order corrections are not small - 
in fact the NLO in a certain region substantially brings down the LO cross section 
and we have explained in the text why this happens. It is clear from our analysis 
that NLO calculations are necessary for UHE neutrino isoscalar scattering processes 
to get a complete picture, not least for isolating the QCD effects from those due to 
possible physics beyond the Standard Model. 
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We have also discussed briefly, observables that are insensitive to higher order 
corrections. Such observables, which are 'stable' under QCD corrections would be 
useful discriminators in the search for new physics. 

At this point, one should also perhaps consider how our results would change 
if one were to use other evolution equations like BFKL [36] or the unified BFKL- 
DGLAP equation of [37]. There have also been attempts to use the BFKL power law 
rise of the structure function and relate this power to the average inelasticity [38]. 
These approaches are relevant in the ultra low x regions that we are discussing. How- 
ever, as has been shown in [34], resumming leading and non- leading ln(l/x) effects 
through the above equations gives neutrino cross sections which are compatible with 
those obtained from the NLO DGLAP framework, thereby reducing potential ambi- 
guities in the extrapolation of the cross section from the (x, Q 2 ) domain of terrestrial 
accelerators (particularly HERA) to the ultra low x region that are explored by UHE 
neutrinos. 

It is possible now to look at the next order in a s (NNLO) to see what effect this 
has at least upto values of x ~ 1CT 4 [40]. This is the subject of a future study. 
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